national  technical 

INFORMATION  SERVICE 


UNITED  ST 

I  r  r  A  1  V 

A  4  ‘  *  ;  A  H  1*1  t  ' 

■  1  < 

NATICK  LAI 

<  n  w  5  I  n  R  1  f  5 

1  r  I 

' 


Best 

Available 

Copy 


Thlt  document  hat  been  approved 
for  public  release  and  sales  Its 
distribution  It  unlimited* 


AO 


TECHNICAL  REPORT 
TR-7 1-32-OSD 


A  GENERAL  THEORY  OP  PARACHUTE  OPENING 


by 


Edward  W.  Ross,  Jr. 


January  1971 


Office  of  the  Scientific  Director 
U.S.  Army  Natick  Laboratories 
Natick,  Massachusetts 


PORE WORD 


The  design  of  parachutes  cannot  yet  be  done  In  a  completely  rational  manner. 
One  of  the  main  reasons  for  this  Is  that  the  opening  process,  during  which 
the  canopy  experiences  Its  most  severe  stresses,  Is  not  well  understood. 

This  report  la  an  attempt  to  clarify  our  comprehension  of  parachute*openlng 
and  systematize  the  estimation  of  Its  main  parameters. 
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ABSTRACT 


A  simple  mathematical  model  of  the  opening  behavior  of  parachutes  is 
presented.  The  model  predicts  the  dreg  and  velocity  as  functions  of  time  and 
also  gives  an  estimate  of  opening  time.  The  model  is  very  general  and  applicable 
to  almost  any  parachute.  The  characteristics  of  a  specific  parachute  enter  the 
model  by  means  of  certain  constants  that  have  to  be  chosen.  Values  of  these 
constants  are  found  for  flat,  circular  parachutes,  partly  by  estimation  and 
partly  by  numerical  experimentation.  The  resulting  model  of  the  opening 
behavior  of  flat,  circular  canopies  gives  reasonable  agreement  with  experi¬ 
ments  for  several  parachutes  of  different  sizes.  Insofar  as  generality  and 
flexibility  are  concerned,  this  model  seems  to  be  superior  to  present  methods 
of  predicting  opening  behavior  of  parachutes. 
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NOMENCLATURE 
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Cgl’  Cg2 

Cml  *  Cm2 
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HD*  HI 
ho»  hl*  h2 
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rn 
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p,  Q 
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*1 

cf 


Maximum  cross-sectional  area  of  an  ellipse 
Drag  Area 

Effective  net  Inlet  area 

Idealised  projected  area 

Semi  axes  of  an  ellipse 

Drag  Coefficient  based  on  frontal  area 

Constants  in  estimates  of  O(V),  (21) 

Constants  In  estimates  of  apparent  mass,  (IV) 
Drag 

Function  defined  In  (7) 

Acceleration,  of  gravity 

Dimensionless  functions  Introduced  in  (10) 
Constants  in  the  formula  (35)  for  Hj(V) 

Total  effective  mass  of  canopy-cargo  system 
Apparent  mass  of  canopy 
Mass  of  fully-opened  canopy 
Mass  of  cargo  plus  parachute 

Number  of  mesh  points  in  numerical  calculations 

Functions  Introduced  in  (8)  and  (13) 

Flat-circular  radius  of  canopy 
2 

u  / 2,  Function  introduced  In  (5) 

Time 

Time  at  which  Inflation  begins 

Time  at  which  full  inflation  Is  reached,  V  - 

Time  at  which  V  -  Vg 
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u  Velocity  of  canopy  and  load 

Up  u»  Velocities  at  start  and  end  of  Inflation 

um  Velocity  of  steady  descant 

V  Volume  of  canopy 

Vj,  Vj  Initial  and  final  volumes  of  canopy 

V8  Vf  -  AV,  (36) 

6}(6|(£|t68  Constants  In  estimate  of  H^,  (33) 

AV  Increment  of  volume  used  In  numerical  Integration 

H  Function  of  ellipsoidal  shape,  (25) 

6  Angle  between  horlsontal  and  parachute  trajectory 

X  Constant  specifying  vent  pull-down,  (31) 

p  Mass  density  of  air 

T  Dimensionless  time  during  steady  descent,  (AO) 

Ts  Dimensionless  time  at  start  of  full  Inflation. 
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1.  Introduction 


In  attempting  to  improve  the  design  of  parachutes,  especially  with  a 
view  toward  launching  from  a  low  altitude,  it  is  desirable  to  understand  and 
predict  the  behavior  of  parachutes  during  opening.  Previous  work  on  opening 
of  parachutes  has  clarified  a  number  of  points,  but  the  analyses  have  been 
either  too  specific  (i.e.  used  experimental  data  for  a  single  canopy),  or  contain 
questionable  assumptions,  or  are  simply  too  cumbersome  for  practical  purposes . 

Since  experimental  results  for  parachutes  are  very  inconsistent,  it  is  point¬ 
less  to  develop  a  highly  refined  theory.  For  the  present  it  appears  that  we 
need  a  relatively  simple  theory  that  gives  roughly  correct  results  for  different 
kinds  of  canopies. 

This  paper  will  describe  such  a  theory.  The  theory  is  formalized  in  a 
Fortran  II  program  whose  demands  on  the  computer  (as  to  speed  and  memory) 
are  quite  modest. 

Hie  procedure  is  to  write  the  governing  equations  for  parachute-opening 
in  a  form  which  is  as  precise  as  possible  without  containing  any  information 
that  is  specific  to  an  individual  parachute.  The  specific  information  is 
contained  in  certain  input -functions  which  occur  in  the  equations,  but  are 
left  in  general  form  at  first.  We  assume  that  these  functions  depend  on 
canopy  volume,  which  is  adopted  as  the  independent  variable  in  place  of  tine. 
General  solutions  in  the  form  of  indefinite  integrals  are  obtained  for  the 
period  of  canopy  inflation.  The  solutions  depend  on  the  input-functions. 

Most  of  these  can  be  estimated  in  advance,  but  the  most  important  one,  which 
prescribes  the  net  flow  rate  into  the  canopy,  is  found  by  trying  various  functions 
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and  choosing  Che  on<>  that  gives  the  beat  agreement  with  experiments.  Por  this 
purpose  the  theoretical  predictions  are  compared  with  testa  on  three  different, 
flat,  circular  canopies. 

2.  Analysis  of  Opening 

During  the  opening  period  we  assume  that  the  direction  of  the  canopy 
axis  Is  approximately  the  same  (Figure  1)  as  that  of  the  canopy  velocity 
vector,  whose  magnitude  Is  called  u.  The  governing  equations  are  taken  to  he 
Conservation  of  Linear  Momentum, 

<Umu)/ch  -  mo3  ©  —  D  (1) 

Conservation  of  Mass. 

dV/dt  *  AjU  (J) 

and  the  Fluid-Drag  Hypothesis 

D  -  U  Aq  (3) 

The  various  quantities  are  defined  In  the  Nomenclature.  The  functions  m,  the 
total  mass  of  the  system,  Aj,  the  effective  net  Inlet  area  of  the  canopy,  Ap, 
the  drag  area,  and  ©,  the  angle  between  horizontal  and  the  velocity  vector, 
contain  the  specific  Information  for  each  parachute.  Ultimately,  we  shall  have 
to  estimate  each  of  these,  but  for  the  present  we  leave  them  In  general  form. 

In  adopting  (1),  (2)  and  (3)  as  our  basic  equations  we  are  neglecting 
a  number  of  factors,  such  as  angular  momentum  and  centrifugal  force  due  to 
curvature  of  the  trajectory.  These  factors  are  extraneous  to  those  that  are 
most  typical  of  the  opening  period,  and  we  wish  to  lfave  them  aside,  although 
w o  recognize  that  their  effects  are  not  always  negligible  In  practice. 
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The  system  (1),  (2)  and  (3)  can  be  reduced  to  a  single  equation  by  Inserting 
(3)  In  (1),  then  dividing  (1)  by  (2)  and  adopting  the  canopy  volume,  V,  as 
independent  variable. 

udCM/dV  =  (*M  S'"  e/A0 ~(f>  AdAH^A)  (« 

In  this  form  ve  assume  that  ra,  0  ,  Aj  and  Ap  are  functions  of  V. 

If  ve  set 

$-  jy-L ,  Si  -  OiVi 


this  can  be  reduced  to  a  linear,  ordinary  differential  equation  of  first  order, 

U  s/d  V)  +  ^/<W)  +  [f  Ao/ 1 m  Al)]}  5 

-s  Vvtt<j  OX*''' 


Such  an  equation  may  always  be  solved  by  finding  a  homogeneous  solution  and  then 
using  variation-of •parameters  to  find  the  particular  solution.  In  viev  of  the 
initial  conditions  the  result  is 


V-Yl  1 


(6) 


where 

p(  (V)/**(Vi)f  e  P(VN)  (T) 

rv 

f’W=fJ  {A0(N')/I>(V')AiW]}<1Y  (8) 

v-v* 

*/hen  s(V)  has  been  found,  the  other  quantities  of  interest  are  determined  by 
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•> 


(9) 


\/l 

u=  (is) 


dV' 

ulTlA^'l 

Formula*  (6)  to  (9)  provide  a  general  solution  of  the  equations  governing 
parachute  opening,  (1),  (2)  and  (3).  The  only  restriction  on  the  validity  or 
these  equations  is  that  the  various  functions  must  be  such  that  the  integrals 
can  be  evaluated. 

It  is  useful  now  to  introduce  three  new  functions  in  place  of  and  Aj. 
Foremost  among  these  is  the  projected  area,  Ap(V).  Two  other  functions,  li^ 
and  Hj,  are  defined  in  terms  of  AD  and  Aj  as  follows: 


t  =  U  i- 


r 

V. 


A0(V)=C0Af^)H^V) 

AP(M)HI(V'I 


We  must  have  Hj(V)>  0  during  opening  or  we  cannot  use  V  as  the  independent 
variable. 

Then  in  terms  of  these  functions  we  may  write  the  solution  during  inflation 
in  the  following  form: 

V 

r  f  un(v) 

'  0  ).  (11) 

w*.  l 


P(V) 

e. 


FCV)  tf(v)  ££ 
r*W)  MjlN')  ApC V) 


(12) 


(13) 
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S  (V)  =  [S4iQ(V)i/  F(v) 

(14) 

Dm--/,C0H0(\n  A  p(V)s(V) 

(15) 

vj(V)  -  \ 

(16) 

v 

t(VU  f  W'OH 1(V>AP(V) 

(17) 

Equations  (11)  to  (17)  determine  the  solution  to  the  parachute  opening  problem, 
In  the  sense  that  we  can  find  u(V),  D(V)  and  t(V)  and,  therefore,  also  u(t), 
D(t)  and  V(t).  Up  to  this  point  the  model  Is  quite  general.  In  the  next 
Section  we  suggest  estimates  of  the  Input  functions,  A^(V),  m(V),  9(V),  H^(V) 
and  Hj(V)  for  the  case  of  a  solid,  flat  canopy. 


3.  The  Input  Functions  for  a  Solid,  Plat  Canopy 

Here  we  wish  to  discuss  the  hitherto  arbitrary  functions,  m,  sin  0  ,  Ap 
Hr.  and  HJf  which  affect  our  solution. 

The  sinplest  of  these  to  estimate  Is  the  total  mass,  m.  We  can  write  m  as 


rv\(V>-=  vtv  4-pV  V  rv\Jy) 


(18) 


Here  mQ  Is  the  constant  mass  of  the  cargo  and  parachute, Is  the  mass  of  the 
Included  air  and  m.(V)  Is  the  apparent  mass  due  to  the  deceleration.  The  apparent 

3 

|3 

mass  Is  rather  difficult  to  estimate  and  Is  often  thought  to  be  unimportant. 

We  shall  assume  a  general  form  suggested  by  an  earlier^  treatment, 

^avv-c^nv/v,)0™-  '  (l9) 

Here  Cmj  and  are  constants  and  is  the  volume  when  the  canopy  Is  fully 
open.  The  final  formula  for  m  is  therefore 


(20) 
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Usually  we  take  CE1  ■  .375  and  ■  .333. 

The  quantity  0  (V)  is  difficult  to  evaluate  for  the  usual  opening  con¬ 
ditions.  Two  extreme  cases  are  easily  handled,  however.  If  the  canopy  velocity 
vector  ie  horizontal  during  the  opening,  then  sin  9  ■  0  and 

Q-O 

This  is  a  very  great  simplification,  and  is  usually  accurate  when  opening  occurs 
soon  after  launching.  At  the  other  extreme,  if  the  velocity  vector  is  vertical 
during  opening,  as  in  the  case  where  a  free  fall  precedes  the  opening,  then  sin  0=1. 
The  estimates  for  the  intermediate  cases  are  quite  uncertain  because  9  really  depends 
in  part  on  the  time  and  not  merely  (as  we  are  assuming)  on  the  volume.  We  shall 
adopt  as  a  general  form 

(a) 

whence  Cg^  ■  0  implies  horizontal  opening,  and  Cg2  *  0,  Cgi  *  TT/2  imply 
vertical  opening.  Usually  we  shall  take  Gg2  *  1  and  Cgi  *  T T  /2,  so  that  9 
varies  linearly  with  volume,  and  the  canopy  attains  a  vertical  orientation  at 
the  end  of  opening. 

The  relation,  Ap(V),  between  the  projected  area  and  the  volume  is  another 
function  that  has  to  be  estimated.  Since  Ap(V)  itself  doe3  not  occur  in  (6)  to 
(9),  but  is  introduced  in  (10),  we  see  that  in  (11)  to  (lM  Ap(V)  always  occurs 
in  combination  with  Hj(V)  or  HD(V).  It  is  not,  therefore,  necessary  that  the 
estimate  of  Ap(V)  be  extremely  accurate,  for  inaccuracies  in  Ap(V)  can  be 
corrected  by  changes  in  Hj  and  Hp.  What  is  required  is  a  simple,  fairly  good 
approximation  that  matches  the  general  features  of  the  true  area-volume  relation, 
but  need  not  be  accurate  as  to  details. 
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With  this  in  mind  ve  derive  Ap(V)  from  an  idealized  picture  of  the 
opening  canopy  in  which  the  canopy  section  (i.e.  the  section  through  the 
axis)  is  taken  as  a  semi-ellipse  whose  arc-length  is  constant  and  equal  to 
the  flat  circular  diameter,  2RC,  see  Figure  2a.  To  find  the  relation  between 
volume  and  area  for  this  ellipsoid,  ve  write  the  well-known  formulas  for 
volume,  V,  and  cross-sectional  area.  A,  of  the  ellipsoid  with  semi -axes 
a  and  b, 

(22) 

A  =  ttol 

(23) 


Since  the  arc-length  from  vent  to  skirt  is  R0,  ve  have 

R„  =  bEOl) 

^  -  l  —  ( ck  /  b  V 

where  E(^  )  is  the  Complete  Elliptic  Integral  of  the  Second  Kind.  The 
relation  between  A  and  V  is  found  parametrically  from  (22)  to  (25), 


where  'l  X  1,  (a  <-<  b)  corresponds  to  the  early  stages  of  opening,  v-j  —  o 
to  a  hemispherical  shape  and  Y)  <  0  to  the  oblate  (a  >  b)  ellipsoidal  shape  in 
the  final  phase.  The  implied  relation  between  A  and  V  is  shown  in  Figure  3. 

The  parametric  formulation,  (26),  of  the  area-volume  relation  is  analytically 
inconvenient,  so  ve  shall  construct  an  approximation  that  is  both  accurate  ami 
convenient.  To  do  this,  we  first  find  the  maximum  value  of  V  as  given  by  (26), 
which  is 


Vwwkx  “  *  •*»  7 


(27) 
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and  occurs  at 


»i  * -t.ru 


cr 


A  -  t  so  l  t? 


•JL 


(?C) 


Since  dV/dA  ■  0  at  this  point,  i.e.  dA/dV  is  infi.iite,  ve  have  no  hope  of 
finding  a  polynomial  approximation  to  A(V)  that  is  uniformly  good  for 
C>  *^iV**M  .  However,  ve  can  find  a  good,  quadratic  approximation  to  its 
inverse  function,  V(A),  by  choosing  the  coefficients  such  that  V  »  0  when 
A  *  0  and  V  attains  its  maximum,  (27),  at  the  A-value  of  (28).  This  approx¬ 
imate  function  V(A)  is  then  solved  by  the  quadratic  formula,  and  the  resulting 
function  is  taken  as  the  projected  area,  Ap, 

£  7 (29) 
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A  slight  modification  of  this  formula  permits  us  also  to  estimate  roughly 
the  area-volume  relation  for  a  canopy  with  a  pulled-down  vent.  Wc  assume  that 
the  vent  is  pulled  down  to  the  skirt  during  the  whole  opening  period,  and 
that  the  section  through  the  axis  consists  of  two  semi -ellipses  instead  of 
one,  see  Figure  2b.  Proceeding  as  before,  ve  find  instead  of  (26) 


A*  -rro-^Le^)] 


-1 


(30) 


We  see  that,  if  in  (26)  V  is  replaced  by  ((."V /  .ViT  ,  we  obtain  (30). 
Thus,  our  approximate  formula  (29)  can  be  generalized  to 

Ap/R>  un-xv*l{.dcs"( (31) 
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where  "X  ■  1  without  a  pulled-down  wont  and  \  ■  1.7  whan  tha  want  la  pulled 
down  to  tha  skirt.  Tha  maximum  volume  (l.a.  volume  In  tha  full-open  config¬ 
uration)  la  given  toy 


V  ft’  =  .  toii/'x 


(32) 


This  completes  tha  estimation  of  Ap. 

In  estimating  Hp(V)  we  adopt  the  viewpoint  that  Cp  Hp(V)  is  the  affective 
drag  coefficient  (tossed  on  projected  area)  at  each  stage  of  opening.  We  give 
Cp  the  constant  value  that  it  has  for  a  fully  open  canopy  and,  therefore, 
regard  Hp(V)  as  a  multiplier  that  corrects  for  the  difference  between  the  drag 
coefficient  of  a  fully  and  partially  open  canopy.  We  shall  adopt  the  following 
form  for  Hp(V): 

(33) 


For  a  solid,  flat  canopy  the  constants  are  given  the  values  CD  ■  1.5  and 

S,’,>  i  **  l 

which  means  that  the  effective  drag  coefficient,  Hj>,  increases  linearly  from 
.73  at  the  start  of  opening  to  1.5  for  the  fully  open  canopy.  This  estimate 
is  based  on  the  idea  that  in  the  early  stages  the  effective  drag  coefficient 
is  roughly  that  (unity)  due  to  the  impact  pressure  acting  over  the  projected 
area,  tout  reduced  somewhat  (to  .73)  because  of  losses  through  the  vent  and 
fabric . 

The  last  of  our  functions,  Hj(V),  is  the  most  difficult  to  estimate. 
Because  of  (2)  and  (10) 

^(V)-  AjOO/AppJ)  (*) 
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l.e.  Hj(V)  Is  Che  rati-  of  Che  effecClve  neC  inlet  ares  Co  Che  projected  area. 
Thus,  Hj(V)  controls  Che  race  at  which  air  fills  the  canopy  and  so  is  affected 
by  many  factors  such  as  the  flow  field  ahead  of  the  canopy,  lorses  due  to 
fabric  porosity,  flow  out  the  vent,  flapping  of  the  gores,  etc.  It  is  not 
easy  Co  find  a  rational  estimate  for  Hj(V),  but  certain  aspects  of  its  behavior 
are  clear  enough*  First,  if  (as  seems  reasonable)  the  canopy  reaches  a  steady 
state  at  the  end  of  opening,  the  net  flow  rate  into  the  canopy  is  then  zero,  hence 

Aj(V^)  =  0  OV\<i  Uj(Vj)  =■  O 

Second,  the  most  important  single  factor  In  determining  the  general  shape 
of  Hj(V)  is  flow  deflection  around  the  canopy.  In  the  early  phases,  because 
of  the  streamlined  shape  of  the  canopy,  it  seem?  plausible  that  relatively 
little  deflection  takes  place.  Later,  as  the  shape  becomes  blunter,  much 
more  deflection  occurs.  This  leads  us  to  expect  that  the  general  behavior 
of  Hj(V)  is  to  decrease  with  increasing  V,  but  it  Is  difficult  to  be  more 
specific  than  this. 

Because  of  the  uncertainty  in  estimating  Hj(V),  it  is  merely  assumed 
that  Hj(V)  can  be  approximated  by  a  cubic  polynomial  of  the  form 

HjCV)--  4-  iv)  +  5  (35) 

The  values  of  the  constants  hQ,  h^  and  hj  were  left  to  be  determined  by  com¬ 
parison  with  experiments. 

4.  Numerical  Treatment 

Equations  (11)  to  (17)  are  taken  as  the  basis  of  a  Fortran  computer 
program  in  which  all  quantities  are  evaluated  successively  at  V  =  Vj,  Vj  +  &V, 

Vt  +  2£V  ....  where  AV  «  tV*.  “T*) / M, 
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and  N  ■  30  usually.  At  each  step  in  V  the  functions  m  (V),  0(V),  Ap(V), 
hD(V)  and  llj (V )  are  computed  from  (20),  (21),  (31),  L3)  and  (35),  respectively, 
ti.cn  the  integrals  (11),  (13)  and  (IT)  are  evaluated  by  Simpson's  Rule  and 
finally  D(V)  and  U(V)  are  calculated  from  (lU),  (lr>)  and  (l6). 

Since  H j(Vf)  ■  0,  the  integrands  in  (11),  (13)  and  (17)  become  singular 
at  V  ■  Vf  and  the  process  has  to  be  terminated  one  step  earlier,  i.e.  at 

v  ■  \  \  •  &V 

In  particular,  this  means  that  ve  cannot  calculate 

but,  must  take  the  last  obtainable  value, 

VW 

as  our  measure  of  opening  time.  This  somewhat  unsatisfactory  procedure  is 
made  necessary  because  we  are  using  V  as  independent  variable.  At  full 
inflation  and  thereafter  (t  ^  tf )  V  remains  constant  am!  cannot  be  used  as 
the  independent  variable. 

If  ve  wish  to  see  what  happens  to  the  system  after  opening  is  completed, 
ve  may  use  time  as  the  independent  variable  for  t  t*  and  continue  the 
calculation,  jetting  M,  9  and  Ap  equal  to  the  values  at  V  ■  V8.  In  this 
phase  of  the  process,  Equation  (l)  is  easily  solvable  in  closed  form, 

VJ  *  3  ’  V*  A 

lit-  (X- 7^)  )  (38) 


where 


u«  -  (39) 
x  =  -t  i'j  m 


I 
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and  T,  and  u#  are  Che  value*  of  T  and  u  when  V  «  Vs.  The  program  carries 
ouC  this  calculation. 

The  very  early  stage  of  opening  require*  some  special  assumptions. 

During  this  period  air  works  its  way  Into  the  gores  and  opens  them  from  the 
skirt  to  the  vent.  The  area-volume  relation,  (31),  derived  from  the  assumption 
of  an  elliptical  section,  Is  not  accurate  In  this  period,  and  the  time  required 
for  this  Initial  filling  is  very  uncertain.  In  this  treatment  we  by-pass 
these  difficulties  by  starting  the  calculations  at  the  end  of  this  period 
and  making  certain  rough  assumptions  about  the  events  preceding  the  start 
of  the  calculation.  These  are  as  fqllows:  first,  the  speed  remains  constant 
during  this  period!  second,  the  elapsed  time  in  this  period  Is 

\  =  O/i)  R./u.  . , 


and,  third,  the  volume,  Vj,  at  the  end  of  this  period  is  calculated  from  the 
approximate  formula, 


\1=  ->s-i  \  <?. 


CSi  fc. 


(42) 


This  formula  was  derived  by  calculating  the  volume  of  the  gores  when  the  vent 
Is  open,  the  lines  straight  from  the  vent  to  the  confluence  point  and  the  gores 
wide  open.  This  calculation  was  carried  through  for  several  different  para¬ 
chutes  and  the  results  approx lma ted  by  the  formula  (42). 

Several  tests  were  run  In  which  the  effects  of  doubling  or  halving  the 
number  of  mesh  points  were  Investigated.  The  final  choice,  N  =  30,  gives  good 
(but  not  excessive)  accuracy  and  short  running  times. 


5.  Calculation  and  Results 

In  this  section  we  discuss  first  the  determination  of  the  function  for 
solid  flat  canopies  under  standard  conditions,  and  then  describe  the  predictions 
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of  the  theory  when  various  parameter-changes  are  introduced. 

The  function  Hj(V)  is  specified  .  j  the  constantn  hc,  h^  and  h2  in 
( 35 ) .  To  find  these  constants ,  many  trial  cases  were  run  with  various 
values  of  hQ,  hi  and  h2  and  for  various  canopy  sizes  and  launch  conditions, 
tome  guidance  was  obtained  from  tests^  on  flow  fields  around  model  parachutes, 
out  much  exploration  was  still  necessary.  The  criterion  for  an  acceptable 
ret  of  constants  vas  that  it  had  to  give  at  least  rough  agreement  with  the 
following  set  of  experimental  results. 

(i)  For  the  C-9  canopy  with  28  ft.  diameter,  ui  ■  170  ft. /sec., 

m^,  ■  7  lb. -sec/ft.  (i.e.  load  ■  225  lbs.),  the  opening  time**  is  about  one 
second  and  the  maximum  drag1  is  about  1200  lbs. 

(ii)  For  the  G-12  canopy  with  6U  ft.  diameter,  recent  unpublished  tests 
with  U*  *  220  ft. /sec.,  ■  70  lb-sec. /ft.  (i.e.  load  »  2250  lbs.)  gave  a 
maximum  drag  of  about  9000  lbs.  No  estimates  of  opening  time  were  made  from 
these  tests. 

(iii)  For  the  G-ll-A  canopy  with  100  ft.  diameter,  ^  ■  200  ft. /sec., 
mo  ■  120  lb. -sec. /ft.  (i.e.  load  ■  3860  lbs.)  the  opening  time  is  in  the 
range  6-8  seconds  and  the  maximum  drag1  is  8000-9000  lbs.  Some  care  is 
necessary  in  interpreting  these  results  because  the  canopy  includes  a  reefing 
line  which  retards  the  early  phases  of  opening  and  allows  some  deceleration  at 
nearly  constant  volume.  Also,  the  aircraft  speed  is  greater  than  that  at  which 
opening  begins  because  the  canopy  and  load  are  launched  with  a  negative  velocity 
relative  to  the  airplane.  These  results  are  probably  appropriate  for  an  air¬ 
craft  velocity  of  about  250  ft. /sec.  or  150  kts. 
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Eventually  It  vas  found  that  the  values 

h*  *  h,  e  .  s  O 

or 

Mt(^=.3ru-(v/v^l 

•J  (1*2) 

give  the  results  summarized  in  Table  1  and  shown  in  Figures  4,  5  and  6.  The 
results  for  maximum  drag  are  roughly  correct  for  the  three  caropies,  and  the 
opening  times  for  the  C-9  and  G-ll  are  also  in  the  right  range.  Accordingly, 
ve  assume  henceforth  that  (1*2)  defines  the  "standard"  behavior  of  the  function 
Hj  (V),  at  least  for  solid,  flat  canopies. 

If  ve  compare  the  drag  curves  qualitatively  for  the  three  canopies  under 
the  conditions  analyzed,  ve  see  that  the  general  shapes  of  the  C-9  and  G-12 
curves  are  similar,  but  for  the  G-ll  maximum  drag  occurs  distinctly  earlier 
than  for  the  others. 

It  is  profitable  next  to  Inquire  whether  the  drag  and  opening  time  are 
morbidly  sensitive  to  changes  in  Hj(V).  This  is  examined  by  introducing 
variations  of  various  kinds  about  the  standard  function  Hj(V)  for  the  C-9 
(see  Figure  7)  and  seeing  vhat  are  the  resulting  changes  in  drag  and  opening 
time  (Figure  8).  In  Figure  7,  curve  B  is  the  basic  curve,  given  by  (1*2), 
and  EP  and  LP  are  pulses  in  the  curve  B  for  Hj.  The  general  conclusion  is 
that  the  results  are  not  excessively  sensitive  to  changes  in  Hj(V)  as  long 
as  Hj(V)  i  0.  If  Hj(V)s£  0  somewhere,  the  opening  process  almost  stops,  and 
the  results  can  be  radically  different  from  cases  where  Hj(V)  i  0. 

The  specific  trends  in  these  results  are  as  follows:  (a)  Increasing 
Hj(V)  causes  an  increase  in  drag  and  a  decrease  in  opening  time,  (b)  The 
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TABLE  1 


Canopy 

Name 

Diam 

ft. 

mo 

lb.sec^/ft1* 

u(l) 

ft/sec 

Dmax 

lbs. 

*s 
sec . 

sec . 

Opening 

Shock 

Factor 

C-9 

28 

T 

170 

1,U00 

1.15 

.51 

6.2 

G-12 

6k 

70 

220 

9,050 

2.60 

1.28 

1*.0 

G-ll 

100 

120 

220 

8,1*60 

6.59 

2.08 

2.2 

Summary  of  Basic,  Computational  Results  for  Three  Parachutes 


TABLE  2 


Canopy 

Diam 

9(V) 

Dmax 

ts 

Name 

ft. 

rad ians 

lbs. 

secs 

C-9 

28 

0 

1290 

1.25 

C-9 

28 

variable 

lUoo 

1.15 

C-9 

28 

TT/2 

1530 

1.12 

G-ll 

100 

0 

7920 

11.  lU 

G-ll 

100 

variable 

81*60 

6.59 

Computational  Results  showing  the  Influence  of 


Gravity  on  Onening  of  Parachutes 
20 


UxlO 
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Figure  4:  Drag  (D),  Velocity  (u)  and  time  (t)  as  Functions  of  Volume  (V)  during  Opening  of  C-9  Canopy 

with  ui  =  170  ft./sec.  and  system  Weight  =  225  lbs.  at  sea  Level 


UxlO 
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Figure  5:  Drag  (I>) ,  Velocity  (u)  and  Time  (t)  as  Function*;  of  Volume  (V)  during  Opening  of  G-12  Canopy, 

with  uj  =  220  ft./sec.  and  system  Weight  -  2,250  lbs.  at  sea  Level 


DxlO 
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Figure  6:  Drag  (D),  Velocity  (u)  and  time  (t)  as  Functions  of  Volume  (V)  during  Opening  of  C-ll-A  Canopy. 

with  uj  _  220  ft./sec.  and  system  Weight  ~  ?.#60  lbs.  at  sea  Level 
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Figure  7:  Various  Functions  Hj(V) 
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nd  Time  (t)  as  functions  of  V  for 


opporent  moss 
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Figure  9:  The  Effect  of  Apparent  Ma  s s _  Drag  CD),  Mass  CM)  and  Time  Ct)  3*  Functions  cf  V 


time 


el  ‘ect  on  the  drag  of  pulse*  in  Hj  is  greater  when  ih>  pulse  occurs  later 
in  the  opening  (assuming  cli< *  pulse  produces  the  sane  percentage  increase  in 
I  at  both  stages). 

The  effect  of  app'irei.t  '.vis  can  be  studied  by  c- .ip*  ring  the  basic 
results  for  the  C-9  with  theme  for  the  case  where  apparent  mass  is  wholly 
s.’glected,  i.e.  in  (20)  w.  t./ixe  ■  0.  Figure  )  :1  the  comparison. 

There  is  little  effect  on  opting  tine,  but  onHtli..  •  -•i  t  ma:.-;  |  reduces 

about  a  10^  increase  n 

The  influence  on  the  solution  of  changes  in  ti  ,  [  i  !.-ag  area.  Ag)  is 
seen  ir.  Figure  10,  where  .•.••suits  are  shown  for  three  nues  of  Si*,  i.e. 

We  see  that  increasi.if  the  drag  area  .increases  botu  me  maximum  drag  and  opening 
time  and  causes  the  drag  ni.xi.muir.  to  occur  et.rlicr. 

To  obtain  an  estimate  of  the  effect  of  g.n  ity  on  the  solution,  we 
examine  what  happens  to  the  C-9  canopy  in  the  simple,  extreme  cases  where 
6(V)  ■  0  and  0(V)  ■  TV/2.  The  results  are  displayed  in  Table  2.  We  see  that 
the  maximum  drag  is  about  8?  less  in  the  former  case  and  9 %  greater  in  the 
latter  case  than  the  standard  value.  Likewise ,  opening  times  are  only  slightly 
affected.  For  the  G-ll  the  case  0  3  0  produces  only  a  small  decrease  in  drag, 
but  a  long  opening  time  (11.1  seconds)  compared  with  the  standard  value  (6.6  sec.). 

The  initial  velocity,  u(l),  has  a  profound  effect  on  the  maximum  drag  and  a 
smaller,  but  still  substantial  effect  on  opc-r  i  time.  In  Table  3  we  see  that 
for  both  the  C-9  and  G-ll,  the  maximum  drag  is  roughly  proportioned  to  [p(l)J2 
and  the  opening  time  to  &(l)  J  Both  these  results  conform  to  our  intuitive 

expectations. 
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TABLE  3 


Canopy 

name 

Diam 

ft. 

lb.  sec2/!*1* 

u(l ) 
ft/aec 

alt. 

ft. 

Dmax 

lbs. 

• 

m 

m  o 
P  o 

m 

C-9 

28 

n 

1 

250 

0 

2,930 

.82 

.1*2 

C-9 

28 

10 

170 

0 

2,060 

.97 

.59 

C-9 

28 

7 

170 

0 

1,L0C 

1.15 

.LL 

C-9 

28 

7 

170 

15,000 

\  ,L6o 

.91 

.63 

0-11 

100 

100 

300 

0 

15.200 

5.5 

.26 

G-ll 

100 

150 

220 

0 

11,000 

5.8 

.38 

G-ll 

100 

100 

220 

0 

8.L60 

6.6 

.32 

G-ll 

100 

100 

220 

15,000 

9,1*00 

5.0 

.L7 

Computational  Results  showing  the  Influence  of  L.nrch  Velocity, 
System  Mass  and  Altitude  on  Opening  of  Parachutes 


TABLE  L 


mQ 

lb/ sec 2/ft** 

u(l ) 
ft/sec 

Dmix5 

lbs 

lbs 

t(U 

sec 

(e) 

*  8 
sec 

.0165 

50 

3.31 

9.L 

.L2 

.215 

.0165 

70 

5.53 

16.1 

.31 

.1L5 

.0165 

85 

7.68 

21.5 

.26 

.121 

Comparison  of 

Theoretical 

((t))  and 

Experiments  1 

((e)) 

Results 

for  the  Opening  of  Parachute  Models 
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Increases  In  load,  m  .  cause  changes  of  th>*  same  s  >rt  as  Increases  In 
<(1),  l.e.  maximum  drag  is  increased  and  opening  decreased.  These 

"esulCs  are  also  displays  I  i>  Table  3.  In  addition,  it  appears  that  Increases 
(n  mQ  cause  the  maximum  drag  to  occur  later  In  the  opening  sequence. 

The  effect  of  altitude  »n  the  opening  behavior  is  shown  In  Table  3. 

For  the  C.9,  if  we  compare  identical  tests  at  sea  i«vel  and  15,000  ft., 
we  find  the  latter  causes  .'out  a  4*4  Increase  in  dr.,};,  a  ?<>*4  decrease  In 
opening  time  and  a  substantial  shift  (delay)  in  the  occurrence  of  maximum 
drag.  For  the  C.ll  the  results  are  roughly  similar  except  that  the  increase 
In  drag  Is  127.  Instead  of  47.. 

Finally,  It  is  Interesting  to  see  what  happens  when  we  apply  our  theory 
to  mode'  tests.  Here  we  compare  our  results  with  tests  of  Heinrich  and  Noreen*^ 
on  three  foot  diameter  models  of  a  C.9  canopy  In  a  vertical  launch  with  a 
load  of  .53  lbs.,  see  Table  4.  Clearly  there  is  a  large  discrepancy  between 
the  theory  and  the  model  tests.  The  theory  rather  consistently  estimates  DMX 
and  ts  as  351  and  210%  of  their  measured  values,  respectively*. 

6.  Discussion 

The  opening  model  proposed  In  this  paper  gives  fairly  satisfactory  results 
for  drag  and  opening  time  of  three  canopies  that  cover  the  range  of  those 
In  common  use.  Also,  it  predicts  successfully  a  number  of  qualitative  trends 
that  are  observable  in  experiments,  such  as  the  drag  Increases  and  opening 
time  decreases  that  occur  when  the  load,  launch  velocity  or  altitude  are 
Increased. 

We  may  make  several  comments  about  the  specific  results  obtained.  First, 
chr  shapes  of  the  drag  curves  are  not  wholly  satisfactory  because  the  drag 


30 


maxima  occur  somevhat  earlier  than  is  usually  observed  in  tests,  especially 
for  the  C-9.  Second,  because  of  the  use  of  V  as  independent  variable,  it  is 
inconvenient  to  reproduce  in  our  model  the  behavior  at  the  end  of  opening. 

It  has  been  observed^ that  solid,  flat  canopies  usually  "over-inflate," 
by  vhich  is  meant  that  Ap  increases  past  the  equilibrium  value  (given 
approximately  by  (28)),  into  the  region  vhere  V  decreases,  before  attaining 
its  equilibrium  value.  In  the  interests  of  simplicity  this  phenomenon  is 
not  incorporated  into  the  model.  The  resulting  inaccuracy  tends  to  shift 
the  drag  maximum  to  an  earlier  time  and  may  partially  account  for  the 
unsatisfactory  shape  of  the  drag  curves. 

We  now  make  some  general  comments  concerning  the  parachute-opening  problem. 
First,  ve  wish  to  emphasize  the  difficulties  in  the  definition  of  "opening¬ 
time"  or  "filling-time."  There  is  a  large  element  of  uncertainty  in  our 
predictions  of  opening-time  due  to  substituting  t8  for  tf.  In  a  sense  this 
is  a  feature  of  the  model  vhich  accords  veil  vith  reality,  because  much 
uncertainty  also  attends  experimental  estimates  of  tf. 

Second,  the  contrast  betveen  the  satisfactory  results  for  canopies 
vith  diameters  in  the  range  28  to  100  ft.  (i.e.  full-size  canopies)  and 
the  very  unsatisfactory  results  for  model  canopies  requires  some  comment. 

It  is  possible  of  course  that  this  is  a  deficiency  of  the  theory.  Hovever, 
vith  the  chosen  constants  the  theory  gives  results  that  are  satisfactory  for 
canopies  vhere  the  largest  diameter  is  3.6  times  the  smallest.  The  diameter 
ratio  of  the  model  to  full-scale  C-9  is 
28/3  -  9-3 

vhich  is  only  2.6  times  the  scale  range  among  the  full-scale  canopies.  It 
seems  unlikely,  if  the  flov  fields  are  similar,  that  the  theory  vould  be  so 
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poor  ;\r  the  model  cancpici. .  We  infer,  therefore,  l  l  .  i  ibably  the  flow 
field  around  the  model  i.uj..tantially  different  from  that  around  the  full- 
scale  canopy.  This  could  lo  mused  by  differing  e‘'h eft-  of  porosity  or 
tending  stiffness  in  tt.c  no.:el  and  full-scale. 

Indeed,  a  very  interesting,  recent  investigati  ;  i6  of  the  effect  of 
flexibility  on  parachute  model  tests  shows  that  me!.. rate  increases  in  flex¬ 
ibility  of  the  models  cuuson  rlzeaule  increase?  i  i  musiF-i- •? ;  opening  times 
ai  d  decreases  in  drag  forces.  !  ur  example,  the  fle\i  l<:  model ^  gave  opening 
times  of  .35,  .21  and  .lb  seconds  for  snatch  vein  : i  t  Lor.  of  50,  70  and  05  ft./ 
sec.  respectively,  instead  of  the  previous  values  (Indie  b)  .215,  .lb5,  .121. 

The  agreement  with  the  theoretical  opening  times  or  iul>lc  b  is  better  but 
still  rather  poor,  but  Table  l6  shows  that  even  the  flexible  model  is  still 
much  stiffer  than  the  full-scale  28-foot  car.opy.  "his  suggests  that  a  con¬ 
siderable  portion  of  the  discrepancy  between  theory  and  experiment  for  small 
models  may  be  attributable  to  differences  in  flexibility. 

In  closing  we  emphasize  that  despite  its  crudities  the  mathematical 
model  of  parachute  opening  given  here  (i.e.  Equations  (11)  to  (17)  together 
with  (20),  (21),  (31),  (33),  (35)  and  (b2))  provides  a  relatively  simple 
method  of  estimating  the  drag  histories  and  opening  times  for  flat,  circular 
canopies.  The  results  agree  reasonably  well  with  experiments  on  three  canopies, 
and  the  computer  program  that  embodies  the  model  is  both  short  and  convenient. 

We  anticipate  that  with  suitable  changes  in  the  constants  of  (33)  and  (35), 
decent  results  may  also  be  obtained  for  ribbon  canopies  and  canopies  with 
pulled-down  vents.  Although  the  model  must  be  employed  with  caution,  it 
appears  to  constitute  an  improvement  over  present  procedures  insofar  as 
generality  and  flexibility  are  concerned. 
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A  simple  mathematical  model  of  the  opening  behavior  of  parachutes  is 
presented.  The  model  predicts  the  drag  and  velocity  as  functions  of  time  and 
also  gives  an  estimate  of  opening  time.  The  model  is  very  general  and  applicable 
to  almost  any  parachute.  The  characteristics  of  a  specific  parachute  enter  the 
model  by  means  of  certain  constants  that  have  to  be  chosen.  Values  of  these 
constants  are  found  for  flat,  circular  parachutes,  partly  by  estimation  and 
partly  by  numerical  experimentation.  The  resulting  model  of  the  opening 
behavior  of  flat,  circular  canopies  gives  reasonable  agreement  with  experi¬ 
ments  for  several  parachutes  of  different  sizes.  Insofar  as  generality  and 
flexibility  are  concerned,  this  model  seems  to  be  superior  to  present  methods 
of  predicting  opening  behavior  of  parachutes. 
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